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We study the vacuum stability of a minimal Higgs portal model in which the standard model 
(SM) particle spectrum is extended to include one complex scalar field and one Dirac fermion. 

These new fields are singlets under the SM gauge group and are charged under a global ?7(1) 
symmetry. Breaking of this U(l) symmetry results in a massless Goldstone boson, a massive CP- 
even scalar, and splits the Dirac fermion into two new mass-eigenstates, corresponding to Majorana 
fermions. The lightest Majorana fermion (w) is absolutely stable, providing a plausible dark matter 
(DM) candidate. We show that interactions between the Higgs sector and the lightest Majorana 
fermion which are strong enough to yield a thermal relic abundance consistent with observation 
can easily destabilize the electroweak vacuum or drive the theory into a non-perturbative regime 
at an energy scale well below the Planck mass. However, we also demonstrate that there is a 
region of the parameter space which develops a stable vacuum (up to the Planck scale), satisfies 
the relic abundance, and is in agreement with direct DM searches. Such an interesting region of 
the parameter space corresponds to DM masses 350 GeV < rUu, < 1 TeV. The region of interest is 
within reach of second generation DM direct detection experiments. 


I. INTRODUCTION 


The conspicuously well-known accomplishments of the 
SU{3)cxSU{2)lxU{1)y standard model (SM) of strong 
and electroweak forces can be considered as the apotheo¬ 
sis of the gauge symmetry principle to describe parti¬ 
cle interactions. Most spectacularly, the recent discov¬ 
ery of a new boson with scalar quantum numbers 

and couplings compatible with those of a SM Higgs has 
possibly plugged the final remaining experimental hole in 
the SM, cementing the theory further. 

Arguably, the most challenging puzzle in high energy 
physics today is to find out what is the underlying theory 
that completes the SM. The overly conservative approach 
to this dilemma has been to assess the consistency of the 
SM assuming a vast desert between the electroweak scale 
Mew 10® GeV and the Planck mass Mpi ~ 10®® GeV. 
The relevant physics of the desert hypothesis is deter¬ 
mined by running couplings into the ultraviolet (UV) us¬ 
ing renormalization group (RG) equations. The behavior 
of the running couplings depends sensitively on the weak 
scale boundary conditions, among which the mass of the 
Higgs boson is perhaps the most critical. The measured 
Higgs mass rriH = 125.5 ±0.5 GeV ms] corresponds to a 
Higgs quartic coupling A close to zero when renormalized 
at energies above A ~ 10®® GeV. 

Strictly speaking, next-to-leading order (NLO) con¬ 
straints on SM vacuum stability based on two-loop RG 
equations, one-loop threshold corrections at the elec¬ 


troweak scale (possibly improved with two-loop terms in 
the case of pure QCD corrections), and one-loop effec¬ 
tive potential seem to indicate uih saturates the mini¬ 
mum value that ensures a vanishing Higgs quartic cou¬ 
pling around Mpi, see e.g. [IHIII- However, the devil is 
in the details. More recent NNLO analyses [IHHSO] yield 
a very restrictive condition of absolute stability up to the 
Planck scale 


niH > 


129.4 ± 1.4 


^mt/GeV- 173.1 


- 0.5 


0.7 

ois{rnz) — 0.1184 
0.0007 


± 1.0th 


GeV. (1) 


On combining in quadrature the theoretical uncertainty 
with experimental errors on the mass of the top (mt) and 
the strong coupling constant (cts), one obtains rriH > 
129 ± 1.8 GeV. The vacuum stability of the SM up to 
the Planck scale is excluded at 2cr (98% C.L. one sided) 
for rriH < 126 GeV [TBH^ . 

The instability of the SM vacuum does not contradict 
any experimental observation, provided its lifetime r is 
longer than the age of the universe Tp. Since the stabil¬ 
ity condition of the electroweak vacuum is strongly sen¬ 
sitive to new physics, from the phenomenological point 
of view it is clear that beyond SM physics models have 
to pass a sort of “stability test” Indeed, only 

new physics models that reinforce the requirement of a 
stable or metastable (but with t > Tp) electroweak vac¬ 
uum can be accepted as a viable UV completion of the 
SM [21H35]. 
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From a theoretical perspective some modification of 
the Higgs sector has long been expected, as the major 
motivation for physics beyond the SM is aimed at resolv¬ 
ing the huge disparity between the strength of gravity 
and of the SM forces. Even if one abandons this hierarchy 
motivation, which does not conflict with any experimen¬ 
tal measurement, the SM has many other (perhaps more 
basic) shortcomings. Roughly speaking, the SM is inca¬ 
pable of explaining some well established observational 
results. Among the most notable of these are neutrino 
masses, the QCD theta parameter, and the presence of a 
large non-baryonic dark matter (DM) component of the 
energy density in the universe. Interestingly, if the new 
dynamics couples directly to the Higgs sector, this may 
induce deviations from the usual vacuum stability and 
perturbativity bounds of the SM. However, beyond SM 
physics models are usually driven by rather high scale 
dynamics {e.g., the neutrino seesaw and the QCD ax- 
ion) , in which case there will be a negligible effect on the 
running of the couplings. A notable exception to this 
is the weakly interacting massive particle (WIMP) DM, 
whose mass scale is constrained to be low if produced by 
thermal freeze-out [36| . 

The scalar Higgs portal is a compelling model of WIMP 
DM in which a renormalizable coupling with the Higgs 
boson provides the connection between our visible world 
and a dark sector consisting of SU (3)c x SU (2)^ x {7(l)y 
singlet fields [57MD] . This is possible because the Higgs 
bilinear $f<I> is the only dimension-2 operator of the SM 
that is gauge and Lorentz invariant, allowing for an inter¬ 
action term with a complex singlet scalar S of the form 

AV = . (2) 


Given that S develops a vacuum expectation value 
(VEV), the Higgs mixes with the singlet leading to the 
existence of two mass eigenstates {hi and / 12 ), which in 
turn open the portal into a weak scale hidden sector. 
Despite its simplicity, in fact, this model offers a rich 
phenomenology, and it provides a simple and motivated 
paradigm of DM. 

In this paper we carry out a general analysis of vacuum 
stability and perturbativity in the SM augmented by a 
Higgs portal with a minimal weak scale hidden sector. 
The layout is as follows. In Sec. [H] we outline the basic 
setting of the scalar Higgs portal model and discuss gen¬ 
eral aspects of the effective low energy theory resulting 
from a minimal hidden sector. In Sec. IHII we confront 
the model with a variety of experimental data, including 
direct DM searches, heavy meson decays with missing 
energy, the invisible Higgs width, as well as astrophysi- 
cal and cosmological observations. In Sec. |IV| we derive 
the RG equations and in Sec. |V] we present the analy¬ 
sis of vacuum stability. Our conclusions are collected in 
Sec. ED 


II. MINIMAL HIGGS PORTAL MODEL 

A viable DM candidate must be stable, or nearly so. 
Stability results from either an unbroken or mildly bro¬ 
ken symmetry in the Lagrangian. A discrete Z 2 sym¬ 
metry is the simplest available symmetry to guarantee 
absolute stability of the DM particle. Under Z 2 the SM 
particles are even while the DM particle is odd [H]. The 
required symmetry may be simply introduced by hand 
into the SM, or, more naturally, may remain after break¬ 
ing of some global continuous symmetry. For example, a 
concrete realization of such a hidden sector could emerge 
when a global U{1) symmetry is spontaneously broken 
by a scalar field with charge 2 under that symmetry, and 
so a discrete Z 2 symmetry arises automatically in the La¬ 
grangian. After spontaneous symmetry breaking, fields 
with an even (odd) charge under the global 17(1) symme¬ 
try will acquire an even (odd) discrete charge under Z 2 . 
Consequently the lightest particle with odd charge will be 
absolutely stable, and thus a plausible dark matter candi¬ 
date. The simplest approach to realize this scenario is to 
introduce one new complex scalar field S and one Dirac 
fermion field ip into the SM. These new fields are singlets 
under the SM gauge group, and charged under U{l)w 
symmetry, such that U{l)w{'fp) = 1 and [/(l)w(5') = 2. 
Spontaneous breaking of a global continuous symmetry 
generates a massless Goldstone boson and a CP-even 
scalar, and splits the Dirac fermion into two new mass- 
eigenstates, corresponding to Majorana fermions. 

The renormalizable scalar Lagrangian density of the 
set up described above is found to be 

V^S -V, (3) 

where 

V = fi^i<l>^<l‘ + fX2^S^S + Ai($^^f + A2(S^Sf + AV (4) 
is the potential and 

2?^ = 9m - *ff2T“IU; - igyYB^ (5) 

is (in a self-explanatory notation) the covariant deriva¬ 

tive. In the spirit of [42] . we write S in terms of two 
real fields (its massive radial component and a massless 
Goldstone boson). The radial field develops a VEV (r) 
about which the field S is expanded 

S=±((r}+r(a:)) (6) 

The phase of S is adjusted to make (o;(x)} = 0. Next, we 
impose the positivity conditions [JD] 

Ai >0, A 2 > 0, A 1 A 2 > -A 3 . (7) 

If the conditions are satisfied, we can proceed to the 

minimization of 0 as a function of constant VEVs for 
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the two scalar fields. In the unitary gauge the Higgs 
doublet is expanded around the VEV as 




-f ^ ) 

V2 [(^}+ Hx) ) ’ 


( 8 ) 


where {4>) = 246 GeV. 

The physically most interesting solutions to the min¬ 
imization of @ are obtained for ( 0 ) and (r) both non¬ 
vanishing 


, ,\2 _ 

" AiA 2 -iAi 


(9) 


and 


2 -Ai/ri + ^Aa/rf 


( 10 ) 


To compute the scalar masses, we must expand the po¬ 
tential Q around the minima ([^ and (10). We denote 
by hi and /i 2 the scalar fields of definite masses, mh^ and 
respectively. After a bit of algebra, the explicit ex¬ 
pressions for the scalar mass eigenvalues and eigenvectors 
read 


and 


with 


= Ai(()-)^-bA 2 (r) 2 -C, 


= Ai ((/))2 + A 2 ( r )2 + C , 


C = 


( 11 ) 

( 12 ) 


\/(Ai(<^)2-A2(r)2)2 + (A3(<(>)(r))2 (13) 


than the Higgs doublet and one in which is lighter. The 
small 6 limit leads to the usual SM phenomenology with 
an isolated hidden sector. 

For the DM sector we assume at least one Dirac field 

/ 

- 

V 2 

As advanced above, we assign to the hidden fermion a 
charge U{l)w{'tp) = 1, so that the Lagrangian is invari¬ 
ant under the global transformation Assuming the 

transformation is local we express tp as 

^(x) = V’'(cc)e*“("). (18) 

Now, after r achieves a VEV we expand the DM sector 
to obtain 

^ • dip' + ip'^'-f ■ dip'''^ , 

— (^>' + ip'^ip"") - V - 

— -{ip'^ip' — ip'^'^ip''') ■ da, 

— (ip'^^ip' + ip'ip''') ■ (19) 

The diagonalization of the ip' mass matrix generates the 
mass eigenvalues. 


ip'^ip s''—-^ipip''S. (17) 


m± =m,p± f{r), ( 20 ) 

for the two mass eigenstates 

iP_ = ^ {iP'^^ - iP') and ip+ = ^ {ip"' + Ip') . (21) 


and 


hi 

^2 


/ cos 9 
sin 0 



(14) 


In the new basis, the act of charge conjugation on ip± 
yields 


ip±=ip±, 


( 22 ) 


Here, 9 G [— 7 r/ 2 , 7 r/ 2 ] also fullfils 

Xsifp) (r) 


sin 29 = 


x/(Ai(0)2-A2(r)2)2 + (A3(</>)(r))2 


(15) 


Now, it is convenient to invert (11), (12) and © , to ex¬ 
tract the parameters in the Lagrangian in terms of mea¬ 
surable quantities: rrih ^, sin 29. We obtain 


Ai — 


2 2 
mt rat 

(1 — cos 26) + (1 + cos 26), 


4(^)^ 


4(<^)2 


2 2 
mt mt 

A 2 = ^(l-cos 20 ) + ^(l + cos 20 ), (16) 


A 3 = sin 20 


2 ( 0 ) (r) 


Note that there are two distinct regions of the parameter 
space: one in which the hidden scalar singlet is heavier 


which implies that the fields ip± are Majorana fermions. 
The Lagrangian is found to be 

= ^^+7 • dip+ + ^0-7 • dip- 

- ^m+ip+ip+ - ^m-ip-ip-, 

% — — 

~ 77t(^+7V'- + ■ da, 

4(r) 

- ^r{Pp+iP++ Pp_ip_), (23) 

where a' = 2a{r) is the canonically normalized Gold- 
stone boson |42] . We must now put r into its massive 
field representation, for which the interactions of interest 
are 

/sin0 - , 7 , ^ /cos0 

=-^— hi^2{W+ip+ + ip-ip-) -^—h2,i 

X (0+^:+ +0-0-). (24) 
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TABLE I. Definition of most common variables. 


-I- 

5 

r 

OL 

H 

ftl,2 

As 

6 
w 

f 


Higgs doublet 
Complex scalar field 
Neutral component of 
Massive CP-even scalar 
Goldstone boson 
SM Higgs boson 
Scalar mass eigenstates 

Quartic coupling between SM and hidden sector 
Mixing angle between hi and h 2 
Lightest Majorana fermion (WIMP) 


w — r coupling constant - see Eq. (231 - 


This leads to 3-point interactions between the Majorana 
fermions and the Higgs doublet. 

All in all, the Dirac fermion of the hidden sector splits 
into two Majorana mass-eigenstates. The heavier state 
will decay into the lighter one by emitting a Goldstone 
boson. The lighter one, however, is kept stable by the un¬ 
broken reflection symmetry. Hence, we can predict that 
today the universe will contain only one species of Ma¬ 
jorana WIMP, the lighter one w, with mass m-a, equal 
to the smaller of m±. Therefore, the dark sector con¬ 
tains five unknown parameters: Wu,, A 2 , 9, and 

/. To facilitate the calculation of the WIMP relic den¬ 
sity, throughout we impose a supplementary constraint 
relating some of these free parameters: Am/wu, ^ 1, 
where Am = |m_|_ — m_| = 2|/(r)|. (The most common 
variables used in this article are summarized in Table |T]) 

A cautionary note is worth taking on board at this 
juncture. It is well known that the spontaneous break¬ 
ing of a global t/(l) symmetry have several disconnected 
and degenerate vacua (the phase of the vacuum expec¬ 
tation value (0|S'|0) can be different in different regions 
of space, and actually we expect it to be different in ca¬ 
sually disconnected regions), yielding dangerous domain- 
wall structure in the early universe |431144j . In the spirit 
of [33], it may be possible to explicitly break the sym¬ 
metry introducing (possibly small) terms in V, such that 
the domain walls disappear before dominating the matter 
density of the universe, while leaving (pseudo-) Goldstone 
bosons and the same dark matter phenomenology 
For simplicity, we restrict our considerations to the po¬ 
tential in (|^, but generalizations are straightforward. 


III. CONSTRAINTS FROM EXPERIMENT 

The mixing of r with the Higgs doublet (j) can be ana¬ 
lyzed in a two-parameter space characterized by the mass 
of hidden scalar m?i. and the mixing angle 0, where i = 1 
for a light scalar singlet (i.e. = mu) and 1 = 2 for 


^ Other approaches, if exceedingly fine-tuned, may offer alternative 
solutions SSHIZ]. 


a heavy one {i.e. mh^ = rnu). We begin to constrain 
this parameter space by using data from DM searches at 
direct detection experiments. 


A. Constraints from direct DM searches 


The wN cross section for elastic scattering is found to 
be 


4 rn'jmlf fp + fn 

TT {rriw + 2 


(25) 


where N = ^{n-\-p) is an isoscalar nucleon in the renor¬ 
malization group-improved parton model [481149j . The 
effective couplings to protons fp and neutrons /„ are 
given by 


fp,n — ^ 


72 rUg 27-^^^ 


(P.n) 


Gq m. 


p,n 


^ , 72 TO, 

q—c,o,t ^ 


(26) 


where Gq is the WIMP’s effective Fermi coupling for a 
given quark species, 




(27) 


with ijjq the SM quark field of flavor q. The first term 
in (26) reflects scattering with light quarks, whereas the 
second term accounts for interaction with gluons through 
a heavy quark loop. The scalar spin-independent form 
factors, are proportional to the matrix element, 

(qq), of quarks in a nucleon. Herein we take [SD| 


fsu = 0.016(5)(3)(1), = 0.014(5)(1^)(1), 

fs, = 0.029(9)(3)(2), = 0.034(9)(13)(2), 

/I, = 0.043(21), = 0.043(21), (28) 


in good agreement with the scalar strange content of the 
nucleon from lattice QCD calculations [51]. The gluon 

scalar form factor is given by = 1 - Y.u,d,s 

For the case at hand. 


-|- 

J p ' J n 

2m% 


~ 0.29 


G„ 


V2r. 


with 


Gq 72 /A 3 (r) 

TUq 2mlml^ ’ 

yielding [5^ 

_1 ml,m% / 0.29A3(r)/ \ 
TT (m^, -I-TOAr)2 ^ TO^TO^ j 


(29) 


(30) 


(31) 
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niw (TeV) 



see Appendixj^for details. We may re-express this result 
in terms of the mixing angle, 


a^N = ( 0 . 29)2 


1 

47r 


{m^ + TTiNY 


X 



sin2 29 . 




(32) 


For 0 <C 1, the upper limits on the nucleon-wimp cross 
sections derived by the various experiments translate into 
upper limits on the mixing angle 


i«i < 


■ TOAf) {(j)) 


m%myj / 


'-/12 


X Q 2g \/^wn{iTLw) ■ 


(33) 


To determine / we require the w relic density to be 
consistent with h'^V,uM — 0.111(6) [53]. In our study 
we consider the interesting case in which < rUu, 
and hence the instantaneous freeze-out approximation is 
valid [54]. In this region of the parameter space, the 
re’s predominantly annihilate into a pair of hi’s or co¬ 
annihilate with the next-to-lightest Majorana fermion, 
producing a scalar hi and a Goldstone boson. All of 
the final state hi subsequently decays into a'. We note, 
however, that for rn-u, « rnn/^ one expects dominant 
annihilation into fermions. We have found that for the 
considerations in the present work, the effective thermal 
cross section can be safely approximated by [54] 


lim {owwVm) « XX-T 

Am/m„->0 oZirm^ 


(34) 


yielding 


1/4 




\/9{xf) Mpifl DM h^ 


( 35 ) 


where Xf = rriw/Tf, g{xf) is is the number of relativistic 
degrees of freedom at the freeze-out temperature Ty, and 
iTihi/iTiw ^ 0-8 [M]- In general for WIMP DM x/ « 20 — 
25 [SS]. The precise relation between the WIMP mass 
and the required Yukawa coupling to attain the relic den¬ 
sity condition is shown in Fig. [l] We note that the mass 
upper limit, TOu, < 74 TeV, is in agreement with the uni- 
tarity limit > 1.7 x 10“®w^ [mu,/(l TeV)]^ [55] . 

which implies rriw <110 TeV [57]. 

Using (33) we can now translate the 90% confidence 
limit on the spin independent elastic WIMP-nucleon 
cross section as obtained by direct detection experiments 
into an upper limit of |0|. In Fig. we show constraints 
on this parameter space from direct dark matter searches. 
For > 8 GeV, the most restrictive constraint comes 
from the LUX experiment [55], whereas for rriw < 8 GeV, 
the most restrictive upper limit is from the SuperGDMS 
low threshold experiment [59]. It should be noted that 
indirect DM searches {e.g. by detecting neutrinos from 
annihilation of captured low-mass WIMPs in the Sun) 
also constrain the WIMP-nucleon elastic scattering cross 
section. However, these searches are in general model de¬ 
pendent. For example, for 100% annihilation into T<r“ 
pairs, the Super-Kamiokande Gollaboration m has set 
the current best upper limit on a^N for WIMP masses 
below 8 GeV. Because of the assumed dominant decay 
into SM fields, this limit cannot be used to further con¬ 
strain the (0,m/i) parameter space. 


B. Constraints from heavy meson decay 

For myj ^10 GeV, searches for heavy meson decays 
with missing energy provide comparable bounds |61H63j . 
In particular, the upper limit reported by the BaBar Gol¬ 
laboration H(T(15') —7 -I- $rp) < 2 X I0“® [M] yields an 
upper bound for the mixing angle, 9 < 0.27 [65] A 
stronger constraint follows from LEP limits on the pro¬ 
duction of invisibly-decaying Higgs bosons <Jzhl<xzH < 
10“'^ [681472] . which implies 9 < 10“2 [^. More restric¬ 
tive constraints come from searches for the rare flavor¬ 
changing neutral-current decay —>• reported 

by the BaBar [73H75] . CLEG [77], and BELLE [^ col¬ 
laborations, as well as limits on K'^ —>■ + ]prp from the 
E787 [75] and E949 experiments [50H55] . The resulting 
excluded regions of the {\9\^mh) plane from all these ex¬ 
periments are compared in Fig. [^ with those from direct 
DM searches. 


C. Constraints from LHC and SN1987A 

Before proceeding we note that additional constraints 
on the (|0|,m^) parameter space can be obtained from 


^ Comparable bounds are obtained from searches for B(T(3S) —i 
7 + [ 66 ] and B{J/'ip ^ 7 -I- fj,) 1671 . 
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FIG. 2. Excluded regions of the (|0|, rrih) parameter space from interactions involving SM particles in the initial state and the 
CP-even scalar in the final state, as well from DM direct detection experiments. The horizontal bands indicate bounds are 
from heavy meson decays with missing energy (no significant excess of such decays over background has been observed yielding 
bounds on the processes T —> 7 / 1 , B"*" —> K'^h, —>■ as well as from LEP limits on the production of invisibly-decaying 

Higgs bosons azh/(^ZH- The diagonal bands represent bounds from DM direct detection experiments (Super-CDMS and LUX), 
for different values of the WIMP mass. Note that all bounds other than the LEP bound can be smoothly extrapolated to the 
smallest rrih ~ 35 MeV stipulated by cosmology. 


limits on Higgs decay into invisible particles and from 
emission of a'-particle pairs in a post-collapse supernova 
core. However, these are not direct constraints as they 
depend also on the quartic coupling of the hidden scalar. 
In particular, since invisible decays reduce the branch¬ 
ing fraction to the (visible) SM final states, it is to be 
expected that B{H —>• invisible) is strongly constrained. 


Indeed B{H —)■ invisible) is known to be less than about 
19% at 95%CL [55H57] . This implies exclusion contours 
in the (|0|, rrih) plane as a function of the free parameter 
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FIG. 3. Bounds on the {\9\,mh) including invisible Higgs 
decays and a' emission in a post-collapse supernova core for 
different assumptions about the value of the quartic coupling 
^2. 


A 2 given by [SS] 


| 6 I(A 2 )| < 1.27 X 10"^ 







In addition, the emissivity of a' due to nucleon 
bremsstrahlung {NN —>■ NNa'a') cannot exceed the lim¬ 
its imposed by SN1987A observations: ta' < 7.324 x 
10“^^ GeV [SH]. For typical supernova core conditions 
(T = 30 MeV and p = 3 x 10^^ g/cm^) it is easily seen 

that IA 3 I < 0.011 ( 5 oo"MeV )^ 1 we can 

translate this limit into a bound on the mixing angle via 


^3 {r){4>) 
mjj-ml ■ 


(36) 


By use of rrih ~ we can express this bound as 


1^1 < 




7.65 m; 


'h 


\ml 


-ml\ 


GeV^ 


(37) 


In Fig. I^we show the exclusion contours for the A 2 = 1 
and A 2 = 0.05. For smaller values of A 2 , the excluded 
regions of the {\0\,mh) plane are dominated by upper 
limits on S-meson decay into invisibles. All in all, for 
= niH, we can conclude from Figs. and that 
2 X 10“^ is a conservative 90% GL upper limit on the 
mixing angle. 


For ^ m//, (33) can be rewritten as 


/ l^l < {4>) rriH V<^wNimiu) 


m 


N 


0.29 


~ 2.7 X 10^ Vcr„fV(m^) GeV. 


(38) 


Dedicated searches for DM candidates serve as an es¬ 
sential component of the LHG physics programme. The 
typical experimental signature of DM production at the 


LHC consists of an excess of events with a single final- 
state partilce X recoiling against large amounts of miss¬ 
ing transverse momentum or energy. In Run I, the AT¬ 
LAS and CMS collaborations have examined a variety 
of such “mono-A” topologies involving jets of hadrons, 
gauge bosons, top and bottom quarks as well as the Higgs 
boson in the final state. In particular, the CMS Collab¬ 
oration has reported very restrictive bounds on the DM- 
nucleon scattering cross section from searches in events 
containing a jet and an imbalanced transverse momen¬ 
tum [^. However, it is important to stress that the con¬ 
tact operator approximation adopted in [SD] only holds if 
the mediator is heavy and can be integrated out [21]. If 
the mediator is light and contributes to resonant DM pro¬ 
duction (as in the minimal Higgs portal model discussed 
herein), the contact approximation fails and the mono-jet 
bounds do not apply. Future LHC 14 mono-A searches 
will also probe vertex operators for which the mediator 
between dark matter and quarks is heavy [211123] ) and 
therefore cannot constrain the Higgs portal model dis¬ 
cussed in this paper. 


D. Constraints from cosmology 


Cosmological observations further constrain the model. 
The earliest observationally verified landmarks - big 
bang nucleosynthesis (BBN) and the cosmic microwave 
background (CMB) decoupling epoch - have become the 
de facto worldwide standard for probing theoretical sce¬ 
narios beyond the SM containing new light species. It 
is advantageous to normalize the extra contribution to 
the SM energy density to that of an “equivalent” neu¬ 
trino species. The number of “equivalent” light neutrino 
species. 


Aeff = 


Pk - 
Pvl 


(39) 


quantifies the total “dark” relativistic energy density (in¬ 
cluding the three left-handed SM neutrinos) in units of 
the density of a single Weyl neutrino 


77r^ 



(40) 


where is the energy density of photons (which by to¬ 
day have redshifted to become the CMB photons at a 
temperature of about ^ 2.7 K) [23] • 

Recent results reported by the Planck Collabora¬ 
tion [95] have strongly constrained the the presence of an 
excess A Aeff above SM expectation: A™ = 3.046 [22]. 
Specifically, the 68% C.L. constraints on Aeff from Planck 
TT, TE, and EE spectra, when combined with polariza¬ 
tion maps (lowP) and baryon acoustic oscillation (BAO) 
measurements are |95] : 


{ 3.13 ±0.32 PlanckTT ± lowP, 

3.15 ±0.23 PlanckTT ± lowP ± BAO, 

2.99 ±0.20 PlanckTT, TE, EE ± lowP, 

3.04 ± 0.18 PlanckTT, TE, EE ± lowP ± BAO. 
















The joint CMB+BBN predictions on N^s provide com¬ 
parable constraints. The 95% C.L. preferred range on 
A^eff when combining Planck data (TT, TE, EE-|-lowP) 
with the helium abundance estimated in [HZ] is iVeff = 
2.99 ± 0.39, whereas the combination of Planck data 
with the deuterium abundance measured in [98| yields 
A^eff = 2.91 ± 0.37 [HS]. (See also [HH|.) In summary, one 
fully thermalize neutrino, AA^eff — 1 , is excluded at over 
3 ct. Models predicting fractional changes of AAeff ~ 0.39 
are marginally consistent with data, saturating the Icr 
upper limit. Models predicting, AAeff ~ 0.57, are ruled 
out at about 2 cr. 

As noted in [42] the Goldstone boson a' is a natu¬ 
ral candidate for an imposter equivalent neutrino. The 
contribution of a' to N^g is AAeff = pa' jPv Thus, tak¬ 
ing into account the isentropic heating of the rest of the 
plasma between the decoupling temperatures, and 
we obtain 


AAeff 


4 


(41) 


where g{T) is the effective number of interacting (ther¬ 
mally coupled) relativistic degrees of freedom at temper¬ 
ature T; for example, g(T,^®®) = 43/4j^ For the particle 
content of the SM, there is a maximum of g{T^?^) = 
427/4 (with > rrit). This corresponds to a mini¬ 
mum value of AAeff = 0.027, which is consistent with 
cosmological observations. However, a fully thermalized 
a', i.e. Ty®® = T^,®® is excluded at 90% C.L.. Note that if 
a' goes out of thermal equilibrium while the temperature 
is just above the muon mass 

AAeff = (4/7)(43/57)^/3 = 0.39. (42) 


This corresponds to a number of equivalent light neutrino 
species that is consistent at the Icr level with current 
data. 

The a' decouples from the plasma when its mean free 
path becomes greater than the Hubble radius at that 
time. The a' collision rate with any fermion species of 
mass irif at or below T is of order [32] 


r(r) 


A§myT’^ 


(43) 


whereas the expansion rate of the universe is of order 


H{T) 


T2 

Mpi 


(44) 


^ If relativistic particles are present that have decoupled from the 
photons, it is necessary to distinguish between two kinds of g: gp^ 
which is associated with the total energy density, and ps, which 
is associated with the total entropy density. For our calculations 
we use g = gp = gs- 


We equate these two rates to obtain 




idee 


mt m 


1/5 


'^2 


A3 my Mpi 


(45) 


Now, taking nif = T = we obtain 


TUh 


(Ai 


m 


H 


(46) 


Substituting the conservative value A 3 = 5 x 10 ^ in 
(46) we have mt ~ 500 MeV. Note that if the a' goes 


out of equilibrium when the only massive SM particles 
left are pairs, AAeff = 0.57. In such a case the 

value of m/i would have to be less than given by ( |46[ ) 
by a factor between (mg/m^)^/^ and (mg/m^)^/^ [42j . 
This sets a lower limit on the mass of the hidden scalar: 
ruh ~ 35 MeV. 


IV. RG EVOLUTION EQUATIONS 


One-loop corrections to (|^ can be implemented by 
making Ai, A 2 , and A 3 energy dependent quantities. The 
positivity conditions of 0 then must be satisfied at all 
energies. 

A straightforward calculation leads to the RG equa¬ 
tions for the five parameters in the scalar potential 


dpi 

dt 

dpi 

dt 

dXi 

dt 


dX2 

dt 

dt 


Pi 

pI 

IGtt^ 

1 

IbTT^ 


r2 . r,P2 \ 9 2 3 2 


12Ai + 6 V/ + 2^A3 - -gl - - 5 ^ 


Pi 




8A2+4^^A3+4/^ , 

P 2 J 

24 A 3 -I- A 3 — 6Y* + -gl + -gy 


10A^ + A^-_/- + 4A2/q , 


2 9 „2 


-glgl. + 12X,Y/ - QX^gl - 3Aig^ ) , (47) 

1 

( 6 A 1 -I- 4 A 2 -I- 2 A 3 -I- 3Y^^‘ — -gl 

+2/q. 


where t = \iiQ and Yt is the top Yukawa coupling, with 


dt 


Yt 

1675^ 


-Y^ - 853 " 



(48) 


and = 'J2mtl{(\)) (see Appendix|B for details). The 
RG running of the gauge couplings fol ow the standard 
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form 


dg^ 

dt 

dm 

dt 

dgy 

dt 


jL 

Wtt^ 

jL 

Wtt^ 

1 

167r^ 


-11 

-f 

4 

-n„ 


7 

4 


3 


l6 

TT^ ’ 

22 


4 

1 ' 


19 

— 

-f 

-rig 

-|- — 


— 

3 

3 ® 

6 


96 


41 


9y 


(49) 


where rig = 3 is the number of generations Finally, 

the running of / is driven by m 


dt 


(50) 


V. VACUUM STABILITY CONSTRAINTS 


We now proceed to study the vacuum stability of the 
model through numerical integration of Eqs. (47), (48), 
(49) and (501. To ensure perturbativity of / between the 
TeV scale and the Planck scale we find from (511, 


/I it-to) \-^/^ 


(51) 


yielding /q < 0.7. For normalization, we set t = 
ln(Q/125 GeV) and tmax = ln(A/125 GeV). Now, us¬ 
ing the SM relation = —2/i^, with ttih — 125 GeV, 
and setting (0)^ = 246 GeV at the same energy scale 
Q = 125 GeV we fix the initial conditions for the pa¬ 
rameters g, and A. Throughout we take the top Yukawa 
coupling renormalized at the top pole mass [ 102 ] . 


A. Light scalar singlet 


We integrate the RG equations from vrih^ = mn and 
impose the initial conditions for Ai. 2,3 by putting the 
observed values into (161 


{4'sm)'^ 


= {<!>? 


, = ruH . 


(52) 


The other quantities in (16) mhi = mh,0 and (r) remain 
free parameters. It is easily seen through numerical inte¬ 


gration of (47), (48 1 , (49) and (50), that there are stable 


vacua up to the Planck scale. However, for those stable 
vacua, the required values of 9 and rrih are excluded at 
90% C.L. 


As an illustration, we note that there is a stable solu¬ 
tion for (r) = 2.8 GeV and ruh = 0.3 GeV, which cor¬ 
responds 9 = 0.01. As can be seen in Fig. this region 
of the parameter space is excluded at 90% G.L. Actually, 
for rrih = 0.3 GeV, it can be shown that the mixing angle 
is bounded from below: 9 > 0.004. The argument is as 
follows. The Yukawa coupling / of the Majorana fermion 
does not suppress the growth of A 2 , but does exactly the 
opposite. This is due to the smallness of / and therefore 


< IbA^/^ in dX 2 /dt. As a result, we can simply set 
/ = 0. The RG equation of A 2 then implies a constraint 
on its boundary value: X2\Q=mH < 0.2 or it blows up 
before reaching the Planck scale. For X2\Q=mH = 0-2, 
we need Xslq^mH < —0.28 to have Ai always positive. 
We note that a positive A 3 only makes Ai grows slower 
and does not help the situation. A smaller A 2 |Q=mjf only 
slows down the growth of A 3 and does not improve the 
stability. In other words, the maximum of A 3 |Q=mjf is 
—0.28. Moreover, from (16) we see that the mixing angle 
decreases monotonically when either X^lQ^mn (when it 
is negative) or A 2 IQ=mif increases. So we reach a mini¬ 
mum angle when X2\Q=TnH = 0-2 and A 3 |Q=mff = —0.28, 
which gives 9 — 0.004. Such a value is excluded at the 
90% C.L. 

Next, we show that for nih > 0.3 GeV, the required 
mixing angle for a stable vacuum up to the Planck scale 
is 0 > 0.004. To this end, we rewrite (16) as 


X2='^2x^+'^{2-2x^). 

4y2 4 y 2 V ^ 


A 3 = 2x 


ml — m 


H 


‘^{4>)y 


(53) 

(54) 


where we have taken x = sin0 and y = (r). Now, since 
ruh < ran by increasing ruh we decrease \m\ — rn?fj\ and 
therefore from (54) we see that xjy increases. Conse¬ 
quently, the term 


m 

4y2 


^2x^ - ^2x‘^ 


4y2 


(55) 


in (53) increases (because it is proportional to x^jy^') 
and therefore the other term (x\/y^ decreases. In other 
words, we have both x/y and y rising and therefore x{9) 
increases with increasing mh- 

It should be noted that a theoretical lower limit on the 
mass of the hidden scalar can be obtained by generaliz¬ 
ing the Weinberg-Linde |1031ll04j bound (see also |105) L 
Herein instead we have used experimental data to deter¬ 
mine such a lower limit. For mt < 0.3 GeV, the previ¬ 
ously derived lower bound on 9 can be relaxed. However, 
for nrih = 35 MeV, we cannot reduce the mixing angle 
to a level consistent with searches for the rare flavor¬ 
changing neutral-current decay K'^ —^ tt"*' -|- ftp with¬ 
out sacrificing vacuum stability, i.e. A 3 ^ 1 is required 
to obtain 9 < 10“^. Moreover, the upper limit set by 
SN1987A observations excludes values of nih < 35 MeV, 
for A 2 < 0.2. As an illustration, in Fig. we show a 
comparison between the 9 behavior imposed by vacuum 
stability and the upper limit on the mixing angle derived 
from (37), fixing the quartic coupling of the hidden scalar 


to the fiducial value that saturates the condition of vac¬ 
uum stability, i.e. X 2 = 0 . 2 . 

We conclude that, for = uih, there are no stable 
solutions up to the Planck scale in the allowed region of 
the parameter space. 
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FIG. 4. Comparison of vacuum stability requirements in 
the (6,mh) plane (blue curve) with the upper limit set by 
SN1987A observations. 


B. Heavy scalar singlet 


For energies below the mass of the heavier Higgs / 12 , 
the effective theory is (of course) the SM. In the low 
energy regime the Higgs sector is given by 

^SM D _ A($t$)2 ^ (56) 


and the RG equations are those of SM. To obtain the 
matching conditions connecting the two theories, follow¬ 
ing [26] we integrate out the field S to obtain a La- 
grangian of the form (56). Identifying the quadratic and 


quartic terms in the potential yields 


2 2 

M = Ml 


M 2 


A 3 

2 A 2 


(57) 


and 


A = Ai 1 - 


>2 


4Ai A^ 


(58) 


respectively. This is consistent with the continuity of 
{ft>SM) ^ {4>)\ namely 




mO 

A 


Q—m^u 


M? 


M2 -^3 

2 A 2 


(1 - nifc) 


Q=mh2 


or equivalently 


(59) 


with {(p) given by (©• The quartic interaction between 
the heavy scalar singlet and the Higgs doublet provides 
an essential contribution for the stabilization the scalar 
field potential [26] . 

When we refer to the stability of Q at some energy Q 
(with the use of the couplings at that scale), we are as¬ 
suming that the field values are at the scale Q. Note that 
the field values are the only functional arguments when 


talking about a potential like ([^, and therefore the ap¬ 
propriate renormalization scale must also be at that scale. 
For A 3 > 0, the third condition in Q could potentially 
be violated only for field values {(p) around mh 2 , regard¬ 
less of the renormalization scale Q [3^. Consequently, 
the region of instability is found to be: 


(r) < mh 2 
Q_ < {<p) < Q+, 
" A 3 


2 _ ”^^2 ^3 
8 A 1 A 2 


i± Wi- 


4Ai A 2 


(60) 


where Q* is some energy scale whe re th e extra positivity 
condition is violated; see Appendix |C|^| Therefore, Q± ^ 
when the extra positivity condition is saturated, that 
is A 1 A 2 = A 3 / 4 . From ( [ 6 ^ it follows that Q± ~ 
when all the A^ are roughly at the same scale. If one 
of the Ai _2 is near zero, then can be ^ but 

this region of the parameter space is constrained by the 
condition Ai _2 > 0. The stability for field values at mh 2 is 
then determined by the potential with coupling at scale 
m/i 2 (instead of Q). Therefore, for A 3 > 0, we impose the 
extra positivity condition in the vicinity of mh 2 - Even 
though the potential seems unstable at Q ^ it is 
actually stable when all the field values are at the scale 
Q. Note that the potential with Xi{Q) can only be used 
when the functional arguments (field valnes (p), (r)) are 
at the scale Q. On the other hand, the instability region 
for A 3 < 0 is given by 


> 


mh2 
•\/ 2 A 2 ’ 


c- < 7 ^ < c+, 
(r) 


A3 

' 2 A 1 


1 ± wi- 


4Ai A' 


X 2 


(61) 


1 A 2 


Q. 


and hence is given by the ratio of {p) and (r), which can 
be reached even with both (p) and (r) being ^ ; see 

Appendix]^ Therefore, for A 3 < 0, we impose the extra 
positivity condition at all energy scales. Note that the 
asymmetry in A 3 will carry over into an asymmetry in 9. 

To solve the system - (47), (48), (49) and (50) - we 
run the SM couplings from 125 GeV up to the mass 
scale 171^2 and use the matching conditions to determine 
(^Sm), which in turns allows one to solve algebraically for 
nihi ■ In Fig. we compare the region of the parameter 
space which contains stable vacua up to the Planck scale 


Note that | |60[ l is where the potential can become negative. If the 
third condition in 0 is satisfied, Q± will be imaginary, which 
implies that the potential is always positive. So we need to make 
sure the third condition is satisfied Q± ~ that the po¬ 

tential can never be negative. On the other hand, we only need 
to consider the third condition in this range as for other (^}, the 
potential is positive regardless of the value of — A 1 A 2 . 
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FIG. 5. The red area shows the allowed parameter space in 
the mho vs. S plane under the vacuum stability constraint 
of Eq. Q, with A = 10^® GeV. The blue areas indicate the 
regions of the parameter space that are not excluded by direct 
DM searches for /o = 0.4, 0.5, 0.7, from light to dark shading. 
The perturbative upper bound is defined by A; < 27r. 


0 



m/,, (TeV) 


FIG. 6. Gomparison of three solutions of stable vacua, with 
identical initial conditions except for /o = 0.4 (red dashed 
line), /o = 0.5 (green dot-dashed line), and fo = 0.7 (blue 
solid line). 
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FIG. 7. The allowed parameter space in the mho vs. & 
plane under the vacuum stability constraint of Eq. Q, with 
A = 10^^ GeV (magenta), A = 10^® GeV (green), and 
A = 10^® GeV (red). The blue areas indicate the regions 
of the parameter space that are not excluded by direct DM 
searches for /o = 0.4, 0.5, 0.7, from light to dark shading. 
The perturbative upper bound is defined by Ai < 27r. 


ify the region of the parameter space with stable vacua, 
see Fig.[^ In summary, the superposition of the blue and 
red areas in Fig. [^indicates the region of the parameter 
space which develops a stable vacuum, satisfies the relic 
density condition, and is in agreement with direct DM 
searches. The interesting region of the parameter space 
comprises WIMP masses 350 GeV < rriw < 1 TeVQ The 
region of interest is within reach of second generation 
DM direct detection experiments, such as DEAP3600, 
DarkSide G2, XENONnT, and DARWIN [rUSlITUg] . 


(red area) with the allowed (blue) bands from direct DM 
searches. From (35) it is straightforward to see that the 
heaviest WIMP satisfying the relic density constraint, 
rriyj = 70 TeV, is near the unitarity limit [56) . How¬ 
ever, one can immediately recognize in Fig. that such 
a WIMP mass exceeds the perturbativity limit, /o < 0.7. 
The maximum WIMP mass that simultaneously satisfies 


the relic density constraint in (35) and the fo pertur¬ 
bativity limit in (51) is = 1 TeV. This maximum 
mass then determines the range of the darker blue band 
in the horizontal axis (rri/i^) of Fig. The LUX up¬ 
per bound on the WIMP-nucleon cross section for elastic 
scattering |58j via (38) sets an upper limit on the mixing 
angle. The allowed values of 6 determine the range of 
the (blue) bands in the vertical axis of Fig. The dif¬ 
ferent blue bands correspond to three fiducial values of 
the Majorana coupling fo = 0.4, 0.5, 0.7. It is impor¬ 
tant to stress that the fo dependence of the RG running 
can be safely neglected; see Fig. It is also impor¬ 
tant to stress that new physics thresholds, which may 
appear near the Planck mass, does not significantly mod- 


VI. CONCLUSIONS 

We have studied the vacuum stability of a minimal 
Higgs portal model in which the SM particle spectrum is 
extended to include one complex scalar field S and one 
Dirac fermion field ip. These new fields are singlets un¬ 
der the SM gauge group and are charged under a global 
17(1) symmetry: U{l)w{'f’) = 1 and U{l)w{S) = 2. The 
spontaneous breaking of this 17(1) symmetry results in a 
massless Goldstone boson, a massive CP-even scalar, and 
splits the Dirac fermion into two new mass-eigenstates 
'0±i corresponding to Majorana fermions. The symme¬ 
try breaking yields naturally a WIMP candidate. Fields 
with an even (odd) charge under the global 17(1) symme¬ 
try will acquire, after symmetry breaking, an even (odd) 


® Curiously, the ATLAS Collaboration has reported a 3 (T excess of 
Higgs pair production HH 7766 for mh2 ~ 300 GeV doe]. 
See also [T07|. 
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discrete charge under a Zi discrete symmetry. While 
the SM particles are all even under Z 2 , the Majorana 
fermions 'i/'i are odd. In such a set up the lightest par¬ 
ticle with odd charge, ip-, will be absolutely stable, and 
hence a plausible WIMP candidate. 

We have shown that interactions between the ex¬ 
tended Higgs sector and the lightest Majorana fermion 
which are strong enough to yield a thermal relic abun¬ 
dance consistent with observation can easily destabilize 
the electroweak vacuum or drive the theory into a non- 
perturbative regime at an energy scale well below the 
Planck mass. However, we have also unmasked a small 
region of the parameter space which develops a stable 
vacuum (up to the Planck scale), satisfies the relic abun¬ 
dance, and is in agreement with direct DM searches. This 
region comprises WIMP masses 350 GeV < mw I TeV. 
The region of interest is within reach of second generation 
DM direct detection experiments. 

Needless to say, here we have considered a minimal 
model to ensure that bounding the parameter space re¬ 
mains tractable. However, our extension of the dark sec¬ 


tor enlarges the parameter space sufficiently to contain 
stable vacua up to the Planck scale. 
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Appendix A 


To determine the value of Gqlrriq we look back at (24) along with the SM Yukawa interaction term, which involves 


the mixing of both scalar fields, hi and / 12 . For interactions of WIMPs with SM quarks, the relevant terms are 


.if = 


m„ cos 6, T , m-„ sin 0, / sin 6* - / cos 9 


-hl,2'lpq1pq - 


{^) (</>) 

The scattering oi a w particle off a quark then gives 

. frug sin 9 cos 9 


h2,l1pq1pq 


-hl^2'4’-i’- + 


h2,iip-ip-. 


M = i- 


JmqXsir) 

i r\ r\ (JJQ 


Uq{p') Uq{p) 


t — m? t — TO? 

, ni,2 (12,1 


Uq{p')Uq{p)u{k')u[k). 


This leads to the identification of the effective coupling 


2Gq 

V2 


JXsir) _ Gq V2fX3{r) 




2 


u{k') u{k) 


(Al) 


(A2) 


(A3) 


Appendix B 

To establish the one-loop RG equations for the parameters of the scalar potential, we first compute the one-loop 
effective potential and then impose its independence from the renormalisation scale. To one-loop level, the scalar 
potential is given by U = -I- where is the tree-level potential and indicates the one-loop 

correction to it. To compute the latter it is useful to re-write the tree-level potential (|^ in terms of the real scalar 
fields: 

4) = 1 and S = [xi -\- ix2) ■ (HI) 

72 V <^3 + y 72 

The particular combination of fields relevant for the calculation are + + and >c^ = >c\ + hence 

(|^ can be rewritten as 

^ 22,^ 2 2, 4 | 1 \ 4 | 1 \ 22 

— 2^1*^ + -^X^p X . 


U(o)(^,r) 


(B2) 
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In the Landau gauge the one-loop correction to the tree-level potential (| B2[ ) reads: 


lnM^(^V^)_ ' 




(B3) 


where Ci are constants that depend on the renormalisation scheme. For the MS scheme, we have Ci = 5/6 for vectors, 
and Ci = 3/2 for scalars and fermions. Next, we expand (B3) and we just keep the contributions from the scalar fields, 
the top-quark, the gauge bosons, and the Majorana fermions, 

1 


= 


647r2 


( 3 ^? 


3' 

+ Gi 

, Q 2 
In —^ 

3' 

-k Tr (u^ 

\ n 3' 

In —^ - 

1 - 12 7)2 


1 ' 


~ 2 



~ 2 

L 2j 

J 

1. 2j 


-k 3Tr M: 


, 5 

In 


- W" 


In 




6_ 

where (in a self-explanatory notation) the field-dependent squared masses are, 

^) = Ml + ) 


2 

G2{^P, >f) = M2 + ■ 


'Hi(p,K) = 


+ 3Aiy)^ -k 


X^ipx 




2 I A3,„2 ; 


^2 “t 3A2>r -k ~2^ 


T^{^) = - {Yt^f , 




2 2 
5yV 

4 V -525V‘M'" 
W^{k) = ^ {fnf . 


-92gY‘f 
2 2 
giv 


(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 


We define the beta functions Pi (i = 1... 3) for the quartic couplings, the gamma functions f^e scalar 

masses, and the scalar anomalous dimensions according to: dXi/dt = /3i, dyp/dt = d/i^/dt = 7 ^ 2 M 2 j 

dtp^/dt = and dn^ jdt = . We then extract the RG equations for the parameters of the scalar potential 

by forcing the first derivative of the effective potential with respect to the scale t to vanish 


^y(i) = ^(y(o) + Al^(i)) = 0 , 
dt dV ' ' 


(Bll) 


keeping only the one-loop terms. After a bit of algebra (|B11|) leads to the following equations: 

1 


2 2 
MI<M 


‘fX 

4 


Pi + 4Ai7,^ — 


IGtt^ 


T'm+27^-l^(^12Ai+2||A3 

MI^^ r r. 1/^0 3 3 Ml \ 

T'm.+27.-i^(^8A2 + 4^A3 

24Ai -k A 3 - 6Y^^ + -^2 + g 5 v + ^525 v 


P2 + 4A27,t — 


Stt^ 


IOA 2 + A 3 — -/ 




1 


Ps + 2A3(7,^ + 7>f) — (6A1A3 + 4A2A3 + 2A3) 


= 0 , 
= 0 , 
= 0 , 
= 0 , 
= 0 . 


(B12) 


Requiring that each term between squared brackets vanishes, we arrive at the RG equations for the parameters of the 
scalar potential. Namely, substituting the explicit expression of the scalar anomalous dimensions [S] 


"fv = 


1 

IGtt^ 


ow 2 9 2 3 2 

- ^52 - ^gy 


and = 


(B13) 


into (B12) we obtain (47). 
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Appendix C 

To explore the impact of the complex singlet scalar on the stability of the Higgs sector we follow 
a tree level scalar potential of the form 




$1$ _ 




s'^s- 




and consider 


(Cl) 


For A 3 > 0, the third term can only be negative when either one of the factors is negative. The parameter space for 
< ((/))^/2 is, in principle, described by the effective potential of the SM (with one Higgs). So herein we only 
consider S < {rY j^. As argued in [5^, the most dangerous region of the field configuration is given by S' = oj^ In 
this region, we have 


y($, 0 ) = Ai(Q) 




A2(Q) 


2 


rf 


A3(Q) 


(<(>)^ 


(C2) 


The couplings are now replaced by their values at some scale Q. We take (^) and (r) to be the physical VEV 
and only the couplings \i run. This is possible in some renormalization scheme (like taking vacuum expectation 
|$| = ((^), |S| = (r) as one of the renormalization conditions, which is satisfied trivially for this particular form of 
potential). Keeping only terms with (r) (since (r) >> (0)), the condition V = 0 can be rewritten as, 


Ai(Q)|$|' 


A 2 (Q)(r )4 A 3 (Q)(r) 


= 0 . 


(C3) 


Next, we assume that A 2 (Q)(r)^ ~ —M 2 (Q)^ is almost unchan ged under the RG flow and remains (ie. we 

assume that Xi does not run by much). Under this assumption (C3) becomes 


Ai(g)i$|4- 


A3(Q)r 

4A2(Q) 


l$|2 + 


G 2 


16 A 2 (Q) 


= 0 . 


(C4) 


The solution to this equation gives Eq. (60); the first condition comes from < (r)^/2 ~ ('/') = 

V2|$| US]- 

For A 3 < 0, we can consider a field configuration with both |$| ^ Q, |S'| ~ Q much larger than (r). The point is 
that we only need to find a configuration in which the stability is violated. In this case, we must keep only the quartic 
term and the potential becomes 


U($,5) = Ai|$|4 + A2|5|4 + A3|$r|^r 


(C5) 


On the one hand, following [26j we can duplicate the procedure to obtain (62). These conditions can be satisfied and 
therefore the vacuum becomes unstable. On the other hand, we can just consider the eigenvalues of the matrix 


Ai ^As 
H A 3 A 2 


(C 6 ) 


In fact, the second approach also tells us why in the case of A 3 > 0, a potential with the form of (C5) is in fact stable. 
The eigenvector with the negative eigenvalue is given by 


( —Xi + A2 + \/X\ — 2 AiA2 + A| + A| ^ 

V "3 ) 

When A| > 4 A 1 A 2 , the first component is negative. So it requires either or \S\^ to be negative, which is impossible. 
As a result, for A 3 > 0 we need to consider a particular field configuration to study the instability. 


® The instability region is defined by both relations Q— < $ < Q+ 
and A1A2 < (2A3)“^, with the couplings evaluated at the scale 
The second relation is more likely to be satisfied at a high 


energy scale, and therefore |$| = Q+ is the most dangerous 
region of the field configuration to reach the instability region, 
i.e. V{^,S) = 0. 
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We now relate the two functional forms of the Higgs potential. At the classical level Q differs from ( |C1| ) by a 
constant; that is the vacuum energy is shifted. In fact (4| has a negative vacuum energy ^ (r)'‘ (again 

neglecting all {(j>) corrections) and the instability requires the potential to be smaller than this negative vacuum energy. 

At a particul ar sc ale Q, all the couplings Xi in (|^ can be replaced by Xi{Q) and so that Q can be rewritten 

in the form of (Cl I with some {4>{Q)), {r{Q)) as a combination of Xi{Q) and Note that we can still adopt our 

previous arguments to consider only the configuration = 0. In this case, 


y(d>,0) = /.2(Q)|ci,|2 + Ai(Q)|<i>|" 


(C7) 


The condition for stability is saturated when 


/r?(Q)|$|^ + Ai(Q)|$|V-A2((r)) (r)^ = 0 . 


(C8) 


Solving (10) we hav e = |A 3 ((r)) (r)^. Now, assuming that all the do not run too much along the RG 

flow we obtain ( |C4[ ). 

When Ai ^2 remains relative away from zero, (60) remains a reasonable approximation for the scale Q± between 


which (j.e., Q- < •\/2|d>| < Q+) the potential can become negative. Note that a naive argument for instability using 
only the quartic potential (which is usually how we get to A| > 4 A 1 A 2 ) is only valid for A 3 <0. As a result, the 
potential can only become unstable in a some very particular field configuration. In this region, however, the effective 
potential is not valid since the field values are far away from the scale Q. 
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